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Coherent control of OH-free radicals interacting with the spin-triplet state of a DNA molecule 
is investigated. A model Hamiltonian for molecular spin singlet-triplet resonance is developed. 
We illustrate that the spin-triplet state in DNA molecules can be efficiently populated, as the 
spin-injection rate can be tuned to be orders of magnitudes greater than the decay rate due to 
small spin-orbit coupling in organic molecules. Owing to the nano-second life-time of OH free 
radicals, a non-equilibrium free energy barrier induced by the injected spin triplet state that lasts 
approximately longer than one-micro second in room temperature can efficiently block the initial 
Hydrogen abstraction and DNA damage. For a direct demonstration of the spin-blockade effect, a 
molecular simulation based on an ab-initio Car-Parrinello molecular dynamics is deployed. 



Organic molecules including DNA have been proposed 
to be used in electronic circuits as molecular conductors 
and nanowires [T] . They are attractive alternatives to in- 
organic semiconductors for a variety of optoelectronic de- 
vices, e.g. photovoltaics. Their long spin lifetime [2] can 
make them desirable for applications in spintronics [3]- 
The unique biological functionality of DNA-molecules in 
carrying the genetic information can be used for bio- 
spintronics. 

We investigate here the possibility to couple coherent 
control methods [4] with bio-spintronics. OH-free radi- 
cals that are products of chemical reactions or ionizing 
radiation [5] are highly reactive and toxic, resulting in 
DNA-damage (single/double strand break), and the de- 
velopment of genetic aberrations, carcinogenesis and ag- 
ing. They are highly unstable diatomic molecules with 
a life-time that is reported within nano-seconds and are 
charged neutral with a magnetic moment that is pro- 
duced by nine electrons with an unpaired spin in the 
outermost open shell. Their high reactivity is thus at- 
tributed to the pairing of opposite spin electrons in the 
open orbital of the free radical with an electron in an 
organic molecule. Fenton chemistry studies used in nu- 
cleic acid foot-printing show that a free radical reaches a 
DNA-molecule by diffusion and removes a hydrogen ion 
from its sugar moiety (e.g., a deoxyribose) [6-8 . The 
most probable sites of attack are identified as H 5 / , H 5 " , 
and H 4 ' of the deoxyribose residue in the backbone of the 
DNA (the nomenclature is adopted following Ref. [H]). 
Because of their short life-time of 1 ns, OH radicals gen- 
erated within 1 nm from the surface of the DNA molecule 
can remove a hydrogen ion and form a water molecule. 

We show here that by optical pumping the HOMO 
electron of both, the OH-free radicals and the DNA, 
can be efficiently spin-polarized. Paramagnetic state of 
OH-radical can be oriented along the non-equilibrium 
direction of the triplet magnetic moment of the DNA 
molecule. Such photo-excited DNA follow a phosphores- 
cence decay to the ground state with slow transition rate, 



hence providing statistically enough triplet excitons to 
block OH-DNA reactivity. The reported triplet life-time 
of DNA nucleotides ranges 1/us at room temperature to 
Is at 77K |2J. The Hydrogen abstraction occurs in 1 ps, 
several order of magnitudes faster than the DNA triplet 
decay process. The disparity in time scale between the 
chemical reaction and DNA-triplet life-time allow us to 
use a ps Car-Parrinello ab-initio molecular dynamics to 
simulate the immunization of DNA against OH-free rad- 
icals. 

Recent studies in pump and probe femto-second laser 
spectroscopy of nucleic-acids in solutions have provided 
valuable information on the optical excitations of nu- 
cleotides and their decay pathways to the ground state 
through radiative and non-radiative channels |5l412|. 
These studies have identified the spin multiplicities of 
the excited states including their life-times and quantum 
yields. As shown in Fig. [T] the major photo-excitation of 
all bases of nucleic-acids occur close to the wavelength 
A = 260 nm (UV) where a direct transition from sin- 
glet ground state So to singlet excited state 1 tttt* is ob- 
served. Unlike III-V/II-VI semiconductors that strong 
SO-coupling has led to direct observation of the optical 
spin excitations in which a laser generates a population 
of photo-carriers strongly spin polarized along a direc- 
tion which depends on the polarization state of the laser 
field [15] ; in organic molecules no direct optical transi- 
tion to spin triplet state has been reported. However, the 
emission spectrum governed by inter-system crossing and 
SO coupling has revealed the spin triplet states, 3 7T7t*, 
with energy below 1 irir*. Consistent with the emission 
spectrum, our ab-initio calculation shows that the low- 
est energy singlet-triplet energy gap (Ax) of DNA-bases 
and DNA-backbone are close to 2.68 eV (corresponding 
to A = 463 nm) for DNA-basc, and 2.38 eV (A = 521 
nm) for sugar moiety. 

We employ a model Hamiltonian consisting of four 
HOMO-LUMO molecular levels as shown in Figure [T] 
The figure shows schematically the direct photo-induced 



2 



'^r*:|TsL 

A ~ 260 nm 



T, 



T n :\U) + \lt 




T +l :\n 



T„ -10" 6 - 



■ls 



(a) Spin single 
ground state 


(S=0) 


UJM<3 — 

hco 02 


TlCQ Q[ 


HOMO — 





(b)Spin triplet (S=l) 
excited state 




B 



B 



FIG. 1: Schematically shows the typical energy level diagram 
(top) and the direct transition of spin singlet-triplet resonance 
(bottom) in the presence of external magnetic field. 



spin singlet-triplet transitions (S = — > 1). The molec- 
ular levels shown in Figure [I] represent the ground and 
excited states of the DNA-molecule. In the absence of 
an external magnetic field, the electronic ground state 
of the DNA-molecule is S = spin-singlet with double 
occupancy of the molecular orbitals. The two-fold de- 
generacy of the molecular orbitals is lifted by a weak 
magnetic field that couples to the electron spin through 
the Zeeman interaction [15 , Ez = gUsS ■ B. Here Ez 
is the Zeeman energy, g is the electron (/-factor [g « 2), 
[Lb is the Bohr magneton (/xg w 5.8 x 10~ 5 eV/Tesla), 
and B is the strength of the external magnetic field along 
quantization axis, z. 

In the Born-Oppenheimer approximation, the 5 = 1 
spin-flip excitation is constructed by removing an elec- 
tron from an occupied So state and putting into an un- 
occupied state with opposite spin. By denoting quasi- 
particle energy levels (electron dressed by interaction) 
by £ n a = £n + S(n, a), the energy of one exciton is Ax = 
£mt — £ni — (m 1\ n i \V\n \.,m t). Here the indices n 



and m label HOMO and LUMO, and a = +1, -1 is used 
for spin up and down. e n is the nth single particle molec- 
ular energy level, E(n,<r) = £„, eocc . [2(n, n'\V\n', n) - 
(n,n'\V\n,n')\ is the electron self-energy and summa- 
tion has carried out over all occupied orbitals, V is 
the electron-electron repulsive Coulomb interaction, and 
Ax = ^triplet - -^singlet is the energy of an exciton relative 
to the So ground state energy. An intra-molecular transi- 
tion accompanied with electron spin flip, can be achieved 
by application of an inhomogeneous perpendicular mag- 
netic field B± = Bj_(r)(x coswii — j/sinwii), assuming 
oj\ > 0. The quasi-particle energy levels shown in Fig. [T] 
form a 4 x 4 block-diagonal Hamiltonian, H . It can be 
decomposed to two 2x2 Hamiltonians Hi and H2, i.e., 
H = Hi® H 2 . Here Hi = h[ 0) + H^ is the two level 
Hamiltonian with the basis |n 4-), and \m t) that form the 



diagonal matrix h[ ^ — e n \\n \){n \. \ + e m ^\m t)(m| | 
with e na = e n — a-^-B. Here 7 = gun/ft is the electron 
gyromagnetic ratio, and the HOMO-LUMO single parti- 
cle energy gap is hu m = E m -f - e ni = (e m - e n ) ~ K'yB. 
Application of B± (r , t) allows a direct transition from 
spin down to spin up (as shown in Figure [T]). The Hamil- 
tonian that describes the coupling of the molecular or- 
bitals with in-plane magnetic field can be expressed as 
Hi = -~lh-fBl m [ejip{iuJit)\n i)(m f | + exp(-iu)it)\m t 
)(n 4, |]. Here £?™ m = (n\B±(r)\m) is the off-diagonal 
matrix element of inhomogeneous in-plane magnetic field 
that allows mixing of HOMO and LUMO. Finally 



Hi = 



£4 

A* exp(— iuiit) 



A exp(iajii) \ 



where A = -fi'yBf' 1 /2. Similarly 



Ho = 



A 



e n i Aexp(—iu>it) \ 
explicit) e m i J 



(1) 



(2) 



In H2 the quasi particle energy gap is hujQ2 — e m \, ~ 
e n + = (e m — e n ) + TljB. Note that the oscillat- 
ing in-plane magnetic field in Hi provides the spin- 
flip resonance with the transition probability given by 
Pni^mt = (\A/h\ 2 /n 2 R )sm 2 (n R t/2). Here £l R = 
vVoi - ^i) 2 + |A/S| 2 . Note that fl R = A/h is the cor- 
responding Rabi frequency at the resonance w i = k>i 
where the optical HOMO-LUMO transition accompany 
with the spin flip-flop shows the highest probability. In 
the last equality we ignored electron-hole Coulomb at- 
traction. Inclusion of V is expected to shift the pre- 
dicted resonance peak to lj± = Ax/h. In contrast, 
because of negative sign of Ui in i? 2 , the oscillating 
in-plane magnetic field does not make a resonance be- 
tween opposite spins in HOMO and LUMO. Therefore 
the spin flip transition induced by B± occurs only be- 
tween I to 4-} an d \m t}> i n accordance with angular 
momentum conservation law. The transition time for 
spin singlet to triplet of a 7r-pulse can be calculated by 

c^/(Bj). The 



{2tt / g)(mc /e£±)/c where £± 
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FIG. 2: Initial state of deoxyribose molecule in the pres- 
ence of OH free radical. Schematically shown the injection 
of photo-generated electrons in deoxyribose with spin polar- 
ization (shown by arrows) along the direction of circularly 
polarized light and external magnetic field. The net mag- 
netic force between two parallel magnetic moments localized 
in OH and deoxyribose is repulsive. This is similar to two 
separated magnetic moments which interact like Heisenberg 
antiferromagnetic exchange coupling. 



steady state population of spin-triplet state tit can be 
calculated by tit = nor r /(r r + t p ) where uq is the den- 
sity of the DNA-base in solution or on the surface-density 
of DNA-bases on the chromatin in nucleus of the cell ac- 
cessible to the solvent (free radicals). As r r >> t p , the 
optically pumped spin-triplet state quickly saturates. A 
laser field with intensity of lW/mm 2 is technologically 
achievable that gives t p « 10~ 7 s. 

Figure [2] schematically shows optically excited 
electron-hole pair (exciton) in the DNA-backbone in the 
presence of spin-polarized OH free radical. Because the 
absorption of a photon alters the local electronic state 
of a DNA-molecule, we confine our simulation to a par- 
ticular segment, e.g., only a part of DNA where the in- 
jected exciton is presumably localized and the optical 
transition takes place, e.g., close to the site of attack. 
Our DNA-backbone model is similar to the deoxyribose 
residue used in other ab-initio calculation such as Ref. 
[13] where the system of interest consists of a DNA nu- 
cleobase modeled by an amino group attached to the de- 
oxyribose in the presence of the OH-radical as shown in 
Figure [2j We selectively choose the site of OH-radical 
attack in the vicinity of Hs< , H5'/ , and H4/ of deoxyri- 
bose residue in the back-bone of DNA because of the 
accessibility of these sites to the solvent [ME] ■ Following 
Ref. [H] we identify the dehydrogenation of the deoxyri- 
bose as a function of spin multiplicity. The ground and 
excited states of the deoxyribose correspond to spin sin- 



FIG. 3: The state of de-hydrogenated deoxyribose by OH 
free radical at t = 0.15 ps (a few time-steps after initial hy- 
drogen abstraction). In this simulation H 4 / on deoxyribose 
is abstracted. The polarization state of the system is spin 
doublet (S =1/2). 



glet (5 = 0), and spin triplet (S — 1) states. As the life 
time of both spin-polarized OH radical and spin-triplet 
exciton in DNA-backbone is much longer than the sim- 
ulation time, we perform CPMD [H] calculat ion under 
restricted spin configuration. We choose two classes of 
spin-restricted calculations, as the total spin along the 
quantum axis is subjected to the constraints S z = 1/2, 
and 3/2, corresponding to doublet and quartet spin con- 
figurations. In both calculations the initial distance be- 
tween OH-radical and DNA-backbone is considered to be 
approximately 1.5 A. 

The initial and final states of the molecules are shown 
in Figures [2]|4] The final configurations of the molecules 
have been obtained after 0.6 ps where the rearrangement 
of the atomic coordinates have been deduced from a dy- 
namical trajectory calculated by CPMD. According to 
our results, a dehydrogenation of the deoxyribose takes 
place around 0.15 ps for a system with S z = 1/2 (total 
spin-doublet) as shown in Figure § This process leads 
to the formation of a water molecule, and H4/ turns out 
to be the site of Hydrogen abstraction in this simula- 
tion. In contrast, as shown in Figure [4] in the quar- 
tet spin configuration the repulsive exchange interaction 
blocks the exchange of hydrogen and hence the chemi- 
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FIG. 5: Evolution of the separation of oxygen atom in OH 
radical with respect to H 4 < in deoxyribose. The dependence 
on the spin multiplicity is seen. The hydrogen abstraction 
occurs approximately at t = 150 fs in the spin singlet state of 
deoxyribose. The process of initial damage to DNA-backbone 
is blocked in deoxyribose due to excessive spin with direction 
parallel to the spin polarization of OH free radical. 



FIG. 4: The state of radio-resistive deoxyribose at t = 0.6 ps. 
The polarization state of the system is spin quartet (S = 3/2) 
induced by circularly polarized light in the presence of weak 
magnetic field. Due to injected localized and polarized photo- 
electrons, the dehydrogenation of DNA-backbone does not 
occur. 



cal reaction. To systematically check the convergence of 
the results, we increased the size of the DNA-molecule 
by adding phosphates and nucleotides and constructed 
a double strand DNA. We used a quantum mechanical 
(QM) approach to simulate the Hydrogen abstraction in 
a small segment of DNA molecule and molecular mechan- 
ics (MM) approach to simulate the vibrational modes of 
the rest of molecule |19j . and found no influence on the 
spin-blocking effect. The details of this study will be 
presented elsewhere. To estimate the energy needed for 
the polarization of the deoxyribose in the absence of OH- 
radicals, we calculated the energy of the ground and ex- 
cited states of the gas-phase deoxyribose in spin singlet 
and triplet multiplicities with the spin singlet-triplet en- 
ergy gap Ax ~ 2.38 eV. This provides an estimate for 
the frequency of the optical gap shown in Figure [TJ which 
is within the range of the visible spectrum of the electro- 
magnetic waves, A = 521 nm, corresponding to green 
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FIG. 6: Kohn-Sham energy as a function of time and spin- 
multiplicity of deoxyribose. An almost flat Kohn-Sham en- 
ergy of spin triplet state of deoxyribose (top) is indication of 
the repulsive interaction between magnetic moments of OH 
free radical and S = 1 exciton due to spin-blockade effect. In 
contrast, a drop in Kohn-Sham energy of spin singlet state of 
deoxyribose (bottom) indicates that the dehydrogenation of 
H 4 / occurs approximately at t = 150 fs. 



light. To estimate the stored magnetic energy due to 
the optical injection of spin, we calculated the energy of 
the gas-phase deoxyribose in the presence of one OH-frcc 
radical with spin doublet and quartet multiplicities. For 
the molecules shown in Figure [2] we find the energy gap 
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-^doublet ~ 3.05 eV. Here the excessive 
magnetic energy which originated from spin-spin repul- 
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sive interactions (which resemble the anti-ferromagnetic 
exchange interaction in the Heisenberg model) can be 
deduced to be Ax — Ax ~ 0.67 eV. This energy can 
be interpreted as the excessive energy barrier due to the 
alignment of the spins in deoxyribose and OH, and is the 
source of the magnetic repulsive force which makes the 
diffusion of OH toward DNA-backbone less likely. This is 
in agreement with the results obtained from CPMD. The 
distance between the oxygen atom in OH radical and H 4 / 
in deoxyribose as a function of time and the spin mul- 
tiplicity of deoxyribose is shown in Figure [5j A drop in 
Kohn-Sham energy seen in spin singlet state that is the 
indication of deoxyribose dehydrogenation, disappears in 
case of spin triplet state. Accordingly dehydrogenation 
of H4/ occurs approximately at t — 150 fs. The corre- 
sponding Kohn-Sham energies are shown in Figure [6] In 
spin triplet state of deoxyribose the spin-blockade effect 
is strong enough that it repels OH free radical. 

In conclusion, we have examined spin-blockade mech- 
anism to propose coherent control of DNA-OH interac- 
tion, in particular from the initial damage induced by 
OH free radicals. To achieve an efficient direct pump- 
ing of triplet states from the singlet ground state we il- 
lustrated that the photo-excitation pumping rate can be 
tuned to be faster than the reported phosphorescence life- 
time, however, an indirect transition via Sq — > 1 irir* — > 
3 7T7r* is an alternative optical pumping pathway that can 
be used. The model calculation was performed using ab- 
initio Car-Parrinello molecular dynamics of a deoxyribose 
residue attached to amino-group in the presence of OH 
free radicals. We provided a detailed description of an 
experimental set up for spin-blockade detection. 
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